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On the semiclassical limit of 4d spin foam models 

Florian Conradjfl and Laurent Freidell] 

Perimeter Institute for Theoretical Physics, Waterloo, Ontario, Canada 

We study the semiclassical properties of the Riemannian spin foam models with 
Immirzi parameter that are constructed via coherent states. We show that, in the 
semiclassical limit, the quantum spin foam amplitudes of an arbitrary triangulation 
are exponentially suppressed if the face spins do not correspond to a discrete geom- 
etry. When they do arise from a geometry, the amplitudes reduce to the exponential 
of i times the Regge action. Remarkably, the dependence on the Immirzi parameter 
disappears in this limit. 
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Loop quantum gravity (LQG) is an approach to canonical non-perturbative quantum 
gravity, where the first-order (or connection) formulation of gravity plays a central role. 
Spin foam models arise from the attempt to construct a corresponding covariant (or path- 
integral) formulation of quantum gravity. In both the canonical and covariant approach, one 
central open issue is the semiclassical limit — the question whether these theories reduce to 
general relativity in suitable semiclassical and low-energy regimes. This problem has been 
explored by many authors and from various angles: for example, by the use of semiclassical 
states [ll- pl, b y the extraction of propagators from spin foam models by numerical 

simulations [lol . 11] and by symmetry reduction fl3]^[l3]- ^^is stage, however, there is no 
conclusive evidence that LQG or spin foam models in 4 dimensions do have a satisfactory 
low-energy behaviour. More tangible results have been obtained in 3 dimensions, where the 
classical and quantum theory 
coupled to point particles [l5| 



are far simpler [15|, ll6|, ll7|, ll8|]. In this case, spin foams were 



19 , and it was found that the semiclassical limit is related to 



a field theory on non-commutative spacetime [18 



Over the last years most investigations in 4 dimensions were focused on a model that 
was introduced by Barrett & Crane (BC) in 1997 2^- It can be constructed by starting 



from a 4d BF theory and by imposing suitable constraints on the 5-field |2ll . |22|. These 



constraints are called simplicity constraints and should restrict the -B-field such that it 
becomes a wedge product of two tetrad one-forms. This procedure for imposing simplicity 
was subject to various criticisms: it was argued, in particular, that the BC model could not 
have a realistic semiclassical limit, since its degrees of freedom are constrained too strongly. 

More recently, two new techniques for constructing spin foam models were introduced 
that open the way to a resolution of this difficulty: the coherent state method 23|, based on 



integrals over coherent states on the group, and a new way of implementing the simplicity 
constraints [13]. These techniques led to the definition of several new spin foam models: 



firstly, a model by Engle, Pereira & Rovelli (EPR) |24j . |25[ |. and later models by Freidel 
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& Krasnov (FK7) [26|] that incorporate any value of the Immirzi parameter 7 7^ 1 and 



reproduce the EPR model for 7 = [26|, |27|. Engle, Pereira, Livine & Rovelli [28l| also 



studied the inclusion of the Immirzi parameter and proposed models (ELPR7) which differ 
from FK7 for 7 > 1. A detailed comparison of the Riemannian models has been performed 
in [iO]. Lorentzian versions of these models have been constructed as well [2^, [2^, 29 . 

In this paper, we focus our study on the set of Riemannian models FK7. The main reason 
for this is the result of j30|, where we showed that each of these models can be written as 
a path integral with an explicit, discrete and local action. We will use this path integral 
representation to analyze the semiclassical properties of the spin foam models FK7. 

As shown in jsO], all known 4d spin foam models with gauge group SO (4) can be written in 
a unified manner. One first introduces a vertex amplitude Ay{j^, le, k^j) which depends on a 
choice of S0(4) representations for each face / of the spin foam, a choice of SU(2) intertwiners 
le for each edge, and a choice of SU(2) representations fee/ for each "wedge" (i.e. each pair 
(e/)). This vertex amplitude isiust the SO (4) 15j symbol with SO (4) representations 
expanded onto SU(2) ones (see [30| for more details). If one sums these vertex amplitudes 
without any constraint one simply obtains a spin foam representation of SO (4) BF theory. 

The spin foam models for gravity arise from two restrictions: firstly, a restriction on the 
SO (4) spins in terms of the Immirzi parameter 7, namely 

J |l-7l 

which implements part of the simplicity constraints^. This implies that there is only one 
free SU(2) spin per face, denoted by jf. The second restriction pertains to the set of SU(2) 
"wedge" representations that one should sum over. It is expressed by the choice of a non- 



trivial measure -Dj^. The cross simplicity constraints require |26l . |28| | that this measure 
should be peaked around k = — for 7 > 1. In the ELPR7 models, this constraint is 
imposed strongly, a la Barrett-Crane, while in the coherent state construction of FK7 it is 
implemented weakly. 

The partition function of the spin foam models is given by a sum over spin foams that 
reside on the dual of a triangulation A and satisfy the above constraints: 

if f 

where 

le,kgf e ef V 

Here, the amplitude W]^{jf) contains the sum over all intertwiner and wedge labels le and 
kef, and can thus be regarded as an "effective" spin foam amplitude for given spins jf. dj 
denotes the dimension of the spin j representation. 

The main focus of our work is to find the semiclassical asymptotics of this effective 
spin foam amplitude. As we will see in the next section, this amounts to determine the 
behaviour of W]^{jf) for large spins. We will find that, in this limit, the effective amplitude 



^ These models are only defined for 7 7^ 1. We also assume that 7 > 0, since a change of sign 7 —7 is 
equivalent to swapping j+ and . 
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is exponentially suppressed if the spin labelling cannot be interpreted as areas of a discrete 
geometry. When the spins do arise from a discrete geometry, on the other hand, and when 
7 > 0, the effective amplitude W]^{jf) is given by the exponential of i times the Regge action. 
It is remarkable that the dependence on the Immirzi parameter drops out. The corresponding 
analysis for the EPR model yields that the exponent vanishes, i.e. the effective action is zero. 

The paper is organized as follows: in section [Tll we review the path integral representation 
that is used to derive the semiclassical approximation. In section IIIIl we define the notion 
of semiclassical limit that we apply in this paper, and present the main result derived in 
the following sections. Section [IV] states the equations which characterize the dominant 
contributions to the semiclassical limit. In sec. |V] we rewrite these equations and project 
them from SU(2) xSU(2) to S0(4). In section IVll we introduce definitions of co-tetrad, 
tetrad and spin connection on the discrete complex. These are needed in section IVIt where 
we show that the solutions to the equations are given by discrete geometries. Finally, in sec. 
IVIIH we put everything together and state the asymptotic approximation of the effective 
spin foam amplitude W]^{jf). 



II. PATH INTEGRAL REPRESENTATION OF SPIN FOAM MODELS 

In this section, we review the path integral representation for the EPR and FK7 models 
derived in ref. 30] and introduce some notations and definitions for simplicial complexes 



and their duals. 

In the following, A denotes a simplicial complex and A* stands for the associated dual 
cell complex. We assume that A is orientable. We refer to cells of A as vertices p, edges i, 
triangles t, tetrahedra r and 4-simplices a. The 0-, 1- and 2-cells of the dual complex A* 
are called vertices v, edges e and faces / respectively. We will also need a finer complex, 
called iSa, which results from the intersection of the original simplicial complex A with the 
2-skeleton of the dual complex A*. This leads to a subdivision of faces / C A* into so-called 
wedges, and each edge e C A* is spht into two half-edges (see Fig. [T]d). We refer to oriented 
half-edges by giving the corresponding pair (ve) or (ev). When an edge in Sa runs from 
the center of a face / to the edge e C df, it is denoted by the pair (/e). A wedge is either 
labelled by a pair ev or by the pair e/, where / is the face that contains the wedge and e is 
the edge adjacent to the wedge that comes first w.r.t. the direction of the face orientation. 

Given Sa and an orientation of its faces /, we define a discretized path integral that is 
equivalent to the spin foam sum ([3]). The variables are spins jf on faces, SU(2) variables 
Ue and Uef on edges and wedges respectively, and SU(2)xSU(2) variables g^e and hg/ on 
half-edges. The set of (g^e^he/) represents a discrete connection on the complex S^- We 
distinguish two types of connection variables, since there are two kinds of half-edges in S^'- 
half-edges (ev) along the boundary df of a face /, and half-edges (e/) that go from an edge 
e in the boundary df to the center of the face / (see Fig. [2]). Given such a connection, and 
for a wedge orientation [eve'f], we can construct the wedge holonomy Gg/ = {G'^j,G~j^), 
where 

Gef = gevgve'K'fhfe . (4) 

The other set of variables {jf,Ue,nef) represent (pre-) geometrical data^. As we will see 
in more detail later, one can think of unit 4-vector normal to the tetrahedron dual 



^ A truly geometrical interpretation is only valid on-shell, when the closure constraint is imposed. 



Figure 1: (a) Face / of dual complex A*, (b) Subdivision of face / into wedges. The arrows 
indicate starting point and orientation for wedge holonomies. 



to e. The spin jf determines the area of the triangle dual to / and rief represents a vector 
normal to this triangle in the subspace orthogonal to Ue- We use these variables to define 
Lie algebra elements X^^ = {X]^ , X]^) G su(2) ©su(2) associated with wedges of Sa- They 
depend on the value of the Immirzi parameter 7 and are given by 



-7 jfUeUefasnJu^^ 



(5) 



Here, cr, denotes the Pauli matrices. The spins jf are arbitrary non-negative half-integers. 
7"'" and 7~ are the integers with smallest absolute value that satisfy 7"^ > and 



2_ 

7" 



7 + 1 
7-1 



(6) 



That is, if 7 > 1, both 7+ and 7 are positive integers, while for 7 < 1, 7 is negative^. In 
the following, we sometimes use the notation 



■7± 



1 7^ I J/ 



(7) 



In the particular cases 7 = and 7 = cx) (corresponding to the EPR and FK model), one 



recovers the usual simplicity relations [20|, |31|, i-e. j 



The action of the path integral is given by 



SlUf, u,, n^j- g.e, he/) = (^(^y ; Gtj) + Six]-- G:^)) 

ejDe 



(8) 



where 



S{X;G) = 2|X|lntr 



^ 1 + 



X 

Jxi 



G 



(9) 



In the last equality, X = X^ai is a SU(2) Lie algebra element, G an SU(2) group element, 
|Xp = X^Xi and the trace is in the fundamental representation of SU(2). Note that by 



In previous papers 



26 



[sot - a different convention was used, where both 7+ and 7^ are positive. This 
entails minus signs in various formulas, depending on whether 7 > 1 or 7 < 1. With the present convention, 
we no longer need to make this distinction, since the minus signs are absorbed into 7" . 
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definition \X]^\ = jj^. This action is invariant under gauge transformation labelled by 
SU(2)x SU(2) group elements Ag, A/, A^, living at vertices, faces and edges of Sa- 



Sev AegeijA^ ^ , hg/ Xeh^fXt^ 



f ' 



n. 



ef 



Up 



(10) 



In order to evaluate this action for a general group element G = PqI + iPg, where Pq = Pc'^i 
and Pq + I-PgP = 1, it is convenient to decompose G into a part parallel to X = X/|X| and 
a part orthogonal to it: 



G 



PgxX 



(^cosei + isinBx) + i - {Pg ■ X)X^ . (11) 



where {Pq x X)* = e^jkP^X^ and cos 6 = P^/y 1 — \Pg x X\'^. The action is in general 
complex, since 



S{X-G) = |X|ln 1 



PgxX 



+ 2i|X|e. 



(12) 



It is important to note that the real part of this action is always negative; Re(S'(X, G)) < 0. 
It is zero only if X is parallel to Pg or equivalently if the Lie algebra element Xcommute 
with the group element G. In this case the action is purely imaginary and has the "Regge" 
form S{X;G) = 2ijX|e. 

As shown in , the spin foam models FK7 introduced in 26| and described in (j3]) can 
be written as 

if f 

where the effective amplitude Wj^ is obtained by integration over all the variables'^ except 
Jf 



(14) 



ejDe 



V, eDv 



ejDe 



III. SEMICLASSICAL LIMIT 

In this section, we define the notion of semiclassical limit that we investigate in this paper, 
and state our main results. We focus our interest on the effective amplitude Wl{jf), which 
depends only on the scalars jf associated to each face. In order to define a semiclassical 
limit we need to reinstate the /i-dependence and introduce dimensionful quantities. The 
spins jf are then proportional to the physical 

The Immirzi parameter enters in the relationship between the discrete bivector field Xj^'^ 
and the dimensionful simple area bivector field A^j associated with the triangle dual to /: 



^ Note that thanks to the gauge symmetry described in pO| there is no need to integrate over the variables 
Ug. The effective amphtude obtained after integration over all variables except jf and is independent 

of Ue- 
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namely,^ 

{16TThG)X]f = ^A,f + -A,j , (15) 

when 7 > 0. 

The simphcity of the area bivector imphes that \A^f\ = \^ef\ = -^Z' where Aj denotes the 
physical area of the triangle dual to /. The relationship (JTSl) can be written (167('hG)'yX]^ = 
(1 ± 7)^^j, which leads to 

^-(7^+7-).,. (16) 

We implement the semiclassical limit by taking h to zero, while keeping the physical 
dimensionful areas Aj fixed. The previous equation (JT6l) tells us that in this limit the spins 
jf are uniformly rescaled to infinity. Thus, the semiclassical regime is reached by taking the 
limit oo of the amplitude W^iNjf) in f|Ti|) . 

Since the action is linear in jf, this corresponds to a global rescaling of the action by 
N. Hence the limit — > cxd is controlled by the stationary phase points of the exponent: 
the integral localizes as a sum over contributions from stationary phase points. Moreover, 
as we have seen, the action is complex with a negative real part. As a result, stationary 
phase points which do not lie at the maximum Re(S'2) = are exponentially suppressed. 
Altogether this means that the semiclassical limit is controlled by stationary points of Sj^ 
which are also maxima of the real part Re(S'2). A more detailed discussion of the asymptotic 
analysis is given in sec. IVIIII Before stating our main result, we have to recall that in the 
continuum the equivalence between gravity and the constrained BF formulation is only 
established if one imposes a condition of non-degeneracy on the B field^. We therefore need 
to distinguish between non-degenerate and degenerate configurations in our analysis. This 
is achieved by splitting the amplitude (fT4l) into two parts^, 

Wl{jf) = wr'{jf) + W2^jf), (17) 

where W^^^{jf) is defined by the integral (fT4ll subject to the constraint that 

\ejjKLX^{i9ee'>X,,f,f''\ > (18) 

for all pairs of wedges (e/) and (e'/') that share a vertex, but do not share an edge. Here, 
gee'>Xe'f' = Qee' > X^i fi and Qee' = Qevdve' ■ The term W^^{jf) denotes the complementary 
integral consisting of degenerate configurations. 

One of the characteristics of 4d spin foam models is the assignment of spins jj to each 
face / of the dual complex A* and of corresponding areas At{if) to each triangle t of A. 
In contrast, Regge calculus is based on an assignment of a discrete metric to the complex, 
defined by lengths associated with each edge £ C A and subject to triangle inequalities. 
The areas At of triangles t dual to faces / are then determined as a function At{h) of the 



The map between bivectors X^j and Lie algebra elements Xe/ — {X^^ai, X ^ai) of su(2) © su(2) is given 
hy^^e^,,Xi';±XO}. 



^ see 



21 



22| for a more detailed discussion of this point and the potential problems due to degenerate 



configurations in the path integral. 
^ see [32.] for an analysis of stationary points of group integrals representing the 6j- and lOj-symbol using 
a similar splitting. 
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edge lengths le- It is well-known 33|, [SJ, l35|] that for an arbitrary assignment of spins jf, 



there is, in general, no set of lis such that Kjf = At{le). The set of areas At determines 
at least one fiat geometry inside each 4-simplex, but the geometries of tetrahedra generally 
differ, when viewed from different 4-simplices. In the following, we will call an assignment 
of spins jf Regge-like if there is a discrete metric k, i C A, such that At{jf) = At{le)- 

Our principal result is that the set of stationary points of the integral (|T4l) which are 
non-degenerate and have a maximal real part, are Regge-like. Moreover, the on-shell action 
is exactly the Regge action! This result relies on the specific realization of the spin foam 
model in terms of the local action ([8]) which is valid for the FK7 version j26| of the model. 



It does not apply to the ELPR7 construction [28|, which is different from FK7 for 7 > 1 
(see ^] for a comparison). 

More precisely, a configuration {jf, Ue, n^f, gev, ^fe) is a solution of the conditions 

dS dS dS dS 



dUef due dgev dhef 



, ReS = , (19) 



and eq. f|T8l) . if and only if the spins jf, f C A*, are Regge-like. In this case, there exist 
edge lengths li, £ C A, such that jf = (7+ + ■y^)~^Af{li) for / = t*. Moreover, for such a 
solution we have, as long as 7 7^ 0, that 

Sl = J2Mh)e>f{k)^Sn{h), (20) 

/ 

where Qf{le) is the deficit angle associated with the face / and Af{le) is the area in Planck 
units. If 7 = 0, the on-shell action vanishes, i.e. S"^ = 0, in agreement with the fact that 
7 = corresponds to a topological theory classically (see (26|). 

It is important to note that the dependance on 7 has disappeared from the functional form 
of the action. This parallels the behaviour of the continuum theory, where the 7 dependence 
drops out classically, once we solve the torsion equation. It also provides a non-trivial check 
on whether the chosen spin foam model captures the right semiclassical dynamics. The 
dependence on the Immirzi parameter arises only at the quantum level as a quantization 
condition on the area^, similar as in canonical loop quantum gravity. 

These results are derived in section IVIII and IVIIH and imply the following statements on 
the effective amplitude W^^^^jf): as N ^ 00, the amplitude W^^'^ (Njf) is exponentially 
suppressed^, if the spins jf, f C A*, do not arise from a Regge geometry. On the other 
hand, if the j/'s are Regge-like, there is a non-zero function CA(j/), independent of A^, such 
that 

Wr^iNjf) ~ ^ (exp (mSn) + cc) (21) 

as N ^ 00. Here, cc stands for the complex conjugate. The number is the rank of the 
Hessian and given by 

TA = 33E -6V -iF, (22) 
with V, E and F denoting the number of vertices, edges and faces of A*. 



The dimensionful area has to satisfy the condition that AtiSirhG) ^(7^ +7 ) ^ is a half-integer. This 
quantization condition becomes invisible in the semiclassical limit h 00. 
That is, the limit iV ^ 00 of N'^W^^'^ {Nj j) is equal to zero for aU neN. 
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Figure 2: Variation of the group variable hjg on the edge fe in the interior of the face /. 

This shows that, in the semiclassical hmit, the effective amphtude W^^'^{jf) is described 
by an effective action, which is the Regge action. If there are several discrete geometries li, 
i G A, for a given set jj, f C A*, one should sum over them in the asymptotic evaluation 
( I2T|) . In the following sections, we prove the above statements and study in detail the 
non-degenerate solutions to eqns. (fT9l) . 



IV. CLASSICAL EQUATIONS 



We will now derive the explicit form of the equations that follow from the conditions 
5S = and ReS = 0. 



A. Variation on interior and exterior edges 

We first consider the variation of the variable /i^^. Since the edge e'f belongs to two 
wedges, denoted ef and e'f, the variation of the action involves only two terms. If e is 
the edge preceding e' along the orientation of the face /, one has Ge/ = gevSve'i^e'fi^fe and 
Ge'/ = gevgi)'e"he"/hje'- Wc cau writc the variation of the action as (see Fig. [2]) 



SS 



2jftT 



Shf,jh%, 



0,(23) 



where we use the abbreviations 



9ee' — 9ev9ve' ' 



{9. 



evJ 



9ve ' 



(24) 



and Xfj: = Xlf / p^, which is independent of 7. 
To write these equations in a more compact manner, let us define the matrix element 



2(1 + ^^/) 



tr Ul + X^f)Gtf 



(25) 



Since 6hh ^ is in the Lie algebra, we conclude from (l23l) that the traceless part of the 
expression in round brackets has to be zero. Moreover, since tr(YgjG^j) = 2, one simply 
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Figure 3: Variation of the group variable gev on the segment ev between faces. 



?ets 



(26) 



We refer to this equation as the interior closure constraint, since it encodes a relation between 
wedges in the interior of the face /. 

Next, we vary a group variable gev on a half-edge ev. This calculation is slightly more 
involved, since the orientation of different faces has to be taken into account. At the edge e, 
four faces /j, i = 1, . . . , 4, intersect. Let be the set of indices i for which the orientation 
of fi is "ingoing" at the vertex v, i.e. parallel to the orientation of the half-edge (ev). In 
these cases, the wedge ho lo no my has the form Gef^ = gevdve^hej^hf.e. Denote by I~ the 
complementary set for which the holonomy is Ge^/j = ge^vOvehef^hf^g.. Then, variation of 
gives 



± 

ev 



6S = tr 



: (27) 

Again, the traceless part of the quantity in round brackets has to be zero. Therefore, 



efi 



j:jtgtiG^.t^.-i 



0. 



(28) 



This equation relates wedges from different faces, so we call it the exterior closure constraint. 



B. Variation of Ue and n^f and maximality 

For the variation w.r.t. rie/, we use the definition ([5]) of and get 

= [5ne/n;/,X+] , 5X- = [uMefn;lu;\ X~] 



(29) 
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The variational equation for ng/ is therefore given by 



Similarly, by varying one obtains 



u 



-1 



Mp = . 



0. 



(30) 



(31) 



The action being complex, its stationarity is not enough to determine the dominant 
contribution to the semiclassical limit. One also has to demand that the stationary points 
are a maximum of the real part of the action. Since 



Re{Sl) = J2 



1 



(32) 



and G]|^ > 0, this is maximal when 




(33) 



Note that the maximum condition implies that the stationarity equations (130113 ip for rig/ 
and Me are automatically fulfilled. Moreover it is important to notice that this relation leads 
to a drastic simplification of the closure constraints, since it leads to the identity: 



(34) 



V. REWRITING THE EQUATIONS 

A. Parallel transport to vertices 

To analyze the variational equations, it is convenient to make a change of variables. 
The original variables Xj^, Ge/ are based at the edge e, which means that under gauge 
transformation they transform as (Xe/, Ge/) (AeX^^A"^, AeGe/A~^). The new variables 
are based at v and defined by parallel transporting the original variables to the nearest 
vertices of the dual complex A*: 

Since every edge e intersects with two vertices v and f', this leads to a doubhng of the 
number of variables. This is compensated by equations that relate variables at neighbouring 
vertices v and v': i.e. 

^Ifiv') = gv'v^lfiv) {gv'vr^ , Ue{v') = g;;,^Ue{v){g^,X^ . (36) 

In terms of the new variables, the interior closure constraint (1261) becomes 

h%gi {Kliv)G%{v)) {h%gf,)-' = h%,gX{K%iv)G^'fiv)) {h%,g^X' ■ (37) 
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Thus, 

G%{v)t%v) = t%{v)Gf,j{v), (38) 

where the edge e' follows the edge e in the orientation of /. Likewise, after conjugation by 
g^^, the exterior closure constraint takes the form 

If we impose, in addition, the maximality constraint flM|) . the closure constraints simplify 
and we remain with the following set of equations: 



(40) 



ee/(f) is a sign factor which is 1 when / is ingoing at v, i.e. oriented consistently with 
the half edge (ev), and —1 otherwise. These equations are supplemented by the simplicity 
constraints 

Kfi^) = (7+J/Ad(ne/(t;))a3 , --f-jfAd{u,{v)n,f{v))as) , (41) 

where 

nef{v) = g^e1^ef{gte)~y (42) 



B. Projection to SO(4) 

In order to solve these equations explicitly it is convenient to project them from 
SU(2)xSU(2) to S0(4) and work purely in terms of vectorial and S0(4) variables. 
The action has the property that 

Slijf, Ue, n^f, -Gef) = Slijf, Ue, n^f, Gef) + 2z7r(7+ + 7")j/ , (43) 

so the weight exp(S'^) projects down to a function of SO (4) if one restricts to configurations 
for which (7+ + 7~)j/ is an integer. We assume from now on that this is the case. 

The projection to SO (4) means that we work with bivectors X^'^ instead of pairs 
(X+, X~), the relation between the two being 

X± = ie,^%fc±Xo,. (44) 

We also associate a unit vector Ue in to each SU(2) element Ue, defined by the relation 



Ue = U^l + iUla,, Ul = l, 



(45) 
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where cTj are the Pauh matrices. To translate the simphcity constraints (HTj) to so(4), it 
is convenient to introduce a fiducial bivector field which is independent of 7 and which is 
simple unlike X'^ . We denote this bivector field by X^f without any subscript 7 and it is 
defined by 

n+ _ v+ _ ^^-X". . (46) 



x2i 



Due to the simplicity constraint (f¥Tl) . u^X^^ + X^^u^ = and \X^j\ = jf. By using the 
identity 

^ (mX+ + X-m) = ■ U)o + i{^X ■ U),a\ (X ■ U)i = XjjU\ {^X)jj = ^ejjKiX''''. 

(47) 

we then find that 



x]f = ^(7+ + i-)x,f + 1(1-'- 7-)(^^e/), i^X^fiv) ■ UMY = 



(48) 



with X(.f{v) ■Xf.fiy) = 2j|. This equation is the discrete version of the simplicity constraints 
in the continuum. Recalling the definition of 7 in terms of 7^, we can write the relation 
between X'^ and X also as 



which shows that for 7 > X plays the role of the dual of the area bivector: 

1 ^Aef 



X 



7+ + 7" SnhG 



(49) 



(50) 



Together with these simplicity conditions, we want to solve the equations (HOl) . When 
written in terms of the 7-independent, simple bivector Xg/, they take the form 



Gef{v)> Xef{v) = Xef{v) , Xef{v')=gy'y>Xef{v), Ue{v') = Qv'vUeiv) 



Xef(v) = X^,f(v 



/3e 



(51) 



> denotes the action of SO (4) generators on bivectors. This and eq. 
of the equations that we will study now. 



HI) are the final form 



VI. DISCRETE GEOMETRY 



In order to find the general solution, we will assume that the bivectors Xf{v) are non- 
degenerate: that is, 

X,j{v)AX,,^,{v)^Q (52) 
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for any pair of faces /, /' which do not share an edge. It turns out that the solutions 
exist only if the set (jf) f is Regge-like. That is, only if there is a discrete metric on the 
triangulation A for which jf is the area of triangles dual to /. As we will see, the unit 
vectors are, on-shell, the normalized tetrad vectors associated with this metric and the 
connection g^i^ is the discrete spin connection for this tetrad. In order to demonstrate these 
statements, we first need to define all these notions on the discrete complex^". 



A. Co tetrads and tetrads on a simplicial complex 

Definition VI. 1 A co-tetrad E on the simplicial complex A is an assignment of vectors 
Ei{v) G Mf^ to each vertex f C A* and oriented edge i C A, £ G a = v* , where the following 
properties hold: 

(i) E_, = -Ee. 

(a) For any triangle t G a = v* , and edges £i,£2,'^3 C t s.t. dt = £i + £2 + £3, the vectors 
Eiiy) close, i.e. 

E,,{v) + Ee,{v) + Ee-M = . (53) 

(Hi) For every edge e = v'v C A* and for any pair of edges £1 and £2 in the tetrahedron t 
dual to e, we have 

E,,{v') ■ Ee,{v') = E,,{v) ■ EM . (54) 

In other words, a co-tetrad E is an assignment of a closed M^-valued 1-chain E{v) to each 
4-simplex a = v* that fulfills a compatibility criterion. In each 4-simplex a* = v, the 
co-tetrad vectors Ei{v) define a flat Riemannian metric g^ by 

giii^iv) = Ei^{v) ■ Ei^{v) , £1, £2 C (J . (55) 

Condition (iii) requires that, for any pair of 4-simplices a = v* and a' = v'* which share a 
tetrahedron r, the metric induced on r by E{v) and E{v') are the same. Thus, the co-tetrad 
E equips A with the structure of a piecewise fiat Riemannian simplicial complex. 

We call a co-tetrad E non-degenerate if at every vertex v = a* (Z A* and for every 
tetrahedron r C o", the span of the vectors Ei{v), £ C r, is 4-dimensional. 

Proposition VI. 2 Given any non-degenerate co-tetrad E, there is a unique S0(4) con- 
nection Q on A* that satisfies the condition 

EfXv') = Vl^>^Et{v) yv'v = ecA\£ce\ (56) 

We call this connection Q the spin connection associated to E. 

Proof Let U{v) and U{v') denote unit normal vectors to Eg-{v), i = 1,2,3, and Ei-{v'), 
i = 1,2, 3, respectively. Choose these unit normal vectors such that 

sgndet (E,,iv'),Ee,iv'),Ee,iv'),Uiv')) = sgndet (e,^{v), Ee,{v), E,^iv),Uiv)) . (57) 
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For previous definitions in the literature, see e.g. [3a, l37|, \3_ 
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Since the tetrad is non-degenerate, we can find a matrix Q^i^ G GL(4) for wliicli 

n,,,Ee^ (v) = Et, {v') ,2 = 1,2,3, n,,,U{v) = U{v') . (58) 

By condition ( 15^ . tliis matrix must be ortliogonal, i.e. ^2^/^, G 0(4). From eq. ( 1571) . we also 
know tliat det f2 = 1, so i7 G S0(4). Suppose now tliere were two matrices ^1,^2 £ S0(4) 
for which 

Q.Ei^v) = EeXv') , Q2EeAv) = EeXv') , 2 = 1,2,3. (59) 
This would imply that 

Q^'Q,E,Xv) = Ee^iv') , 2 = 1,2,3, (60) 
and hence ^2 = ^i- Therefore, the group element Qjj'v is unique. H 

Together, the closure condition (!53|) and equation (!56|) can be regarded as a discrete analogue 
of the torsion equation DE = 0. 

In analogy to the continuum, we can define the concept of a tetrad. This definition makes 
heavy use of the duality between A and A*. To describe the relation between tetrad and 
co-tetrad, it is, in fact, convenient to formulate everything in terms of the dual complex A*. 
Each 4-simplex a is dual to a vertex f of A*, i.e. a = v*. By deleting a vertex p in a, we 
obtain a tetrahedron r. This tetrahedron r is, in turn, dual to an edge e. Thus, the choice 
of a 4-simplex a and a vertex p G a defines an edge e at the dual vertex v. Conversely, a 
pair (t>, e) can be used to label a vertex p of the triangulation. Two different pairs (f , ei) 
and {v',e[) correspond to the same vertex provided that 1) (f,f') = e is an edge of A* and 
that 2) (ei, e, e[) are consecutive edges in the boundary of a face of A*. 

Since vertices of the triangulation correspond to pairs (f,ei), edges i = [P1P2] C a of A 
correspond to triples {v, Ci, 62). We can use this to translate the notation for the co-tetrad 
to the dual complex: instead of denoting the co-tetrad by Ei{v), we can write it as Ee^e2{v)- 

In this notation, the defining relations for the co-tetrad appear as follows: 

Eee' (V) = -Ee'eiv) , E^.e, (v) + E^.e, (v) + E^.e, (v) = . (61) 

Similarly, the equation for the spin connection becomes 

n,E,,,,{v) = E,,,,{v') , (62) 

where e = {vv') is an edge of A* and (ej,e^)j=i^2 are pairs of edges such that (ej,e, are 
consecutive edges in the boundary of a face. Note that there are always four such pairs for 
a given edge e. 

When stating relations between co-tetrad and tetrad, it is also convenient to define 
an orientation for each 4-simplex. By definition, a local orientation of A is a choice of 
Z2-ordering of vertices for each 4-simplex a. Such an ordering is represented by tuples 
[Pi; ■ ■ ■ 5^5]- Two Z2-orderings that differ by an even permutation are by definition equiv- 
alent. Two Z2-orderings that differ by an odd permutation are said to be opposite and 
we write [pi,P2 - ■ ■ ,P5] = — [P2,Pi ■ ■ ■ ,^5]- By duality it is clear that a local orientation is 
equivalent to a choice of Z2-ordering [ei, ■ ■ ■ ,65] of edges of A* meeting at v. With this 
orientation we can also define a correspondence between edges ei and oriented tetrahedra 

Given a choice of local orientation of A, one says that two neighboring 4-simplices a, 
a' that share a tetrahedron r are consistently oriented, if the orientation of r induced 
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from (J is opposite to the one induced from a'. Namely, if o" = [poiPi-, - ■ ■Pi\ and a' = 
— [Pq,Pi, ■ ■ -Pi], they induce opposite orientations on the common tetrahedron r = [pi, ■ ■ - pi] 
and are therefore consistently oriented. The triangulation A is said to be orientable when 
it is possible to choose the local orientations such that they are consistent for every pair 
of neighboring 4-simplices. Such a choice of consistent local orientations is called a global 
orientation. 

From now on and in the rest of the paper we assume that we work with an orientable 
triangulation and with a fixed global orientation. 

Definition VI. 3 For a given non-degenerate co-tetrad E on A, the associated tetrad U is 
an assignment of vectors Ue{v) G M'^ to each vertex v and (unoriented) edge eDv such that 

Uel{v)El„^,{v) = 6e"e " ^e'e (63) 

for all e',e" D v. 

These conditions specify the tetrad U uniquely, as we show in appendix|Al The orthogonality 
relation (163|) is the discrete counterpart of the equation E^^iEj = in the continuum. 
Based on (!63l) . we can derive a number of useful identities satisfied by a tetrad: 

Proposition VI. 4 At any vertex v C A*, the tetrad vectors Ue{v) close, i.e. 

J2u,iv) = 0. (64) 

For a tuple [ci . . . 65] of edges at v, we can express the discrete tetrad explicitly in terms of 
the discrete co-tetrad and vice versa: 

UeM = ^ {E,,M A E,,,,{v) A E,,,,{v)) (65) 

and 

Ee.eM = V,{v) * iUe,{v) A A Ue,{v)) , (66) 

where V4(f)/4! is the oriented volume of the 4-simplex spanned by the co-tetrad vectors: 

V^{v) = det{E,,,,{v),...,E,^,,{v)) . (67) 

For bivectors, one has the relation 

^ iE,,,,{v) A E,.,,,iv)) = V,iv) iUe,iv) A U^^iv)) ■ (68) 

The norm of Ue is proportional to the volume V3(e)/3! of the tetrahedron orthogonal to Ue'- 

\UAv)\ ^ (69) 

The determinant of the tetrad vectors equals the inverse ofV^^iv): 

^ dei{U,,{v),...,U,,{v)). (70) 



In this proposition we have set [*(-Ei A ■ ■ ■ A -En)]/i.../4_„ = ei^...i^E[^ " ■ ■ ■ E^^ . These state- 
ments are proven in appendix lAl 
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VII. SOLUTIONS 



With the help of the previous definitions, we will now determine the solutions to the 
equations (ISTI) . the simplicity constraints ( l48|) and the non-degeneracy condition ( l52l) . For 
the proofs it is practical to denote the oriented wedge (e/) by an ordered pair of edges (ee') 
which meet at v. The order (ee') refers to the fact that e and e' are consecutive w.r.t. the 
orientation of the face. Note that the interior closure constraints Xf,f{v) = Xe'f{v) mean 
that there is only one bivector per face / and vertex v. Hence we can denote the bivectors 

by Xee'{v) = Xef{v) = Xe'f{v). 

Proposition VII. 1 Let {j f,nef,Ue, gev, hef) be a configuration that solves eqns. ( f37]) . with 
the bivectors defined by the simplicity condition [4^ - Then, there exists a co-tetrad E such 
that for any vertex v and tuple [ei . . . e^] of edges at v 

Xe,es{v) = € ^ {Ee.eM A E^^eM) ■ (71) 



The factor e is a global sign, and the ^-volume V^iiv) is given by equation ( 61). This co-tetrad 
is unique up to inversions Ei{v) — > —E^{v), i C v* . 

Equivalently, the bivectors can be expressed by the associated tetrad U , namely, 

X,,,{v) = eV^{v) {Ue{v) A UAv)) (72) 

for any pair of edges e,e' D v. 

Given this co-tetrad E and tetrad U, the variables {jf,nef,Ue, gev) o,re determined as 
follows: the spin jf is equal to the norm of the bivector -kXef{v), and hence Regge-like. The 
group elements g^'v = gv'egev are, up to signs ee, equal to the spin connection for the co-tetrad 
E, i.e. 

gv'v = ^e^v'vi = ±1 • (73) 

For a given choice of the holonomy g^v on the half-edge ev, the group element Ue is deter- 
mined, up to sign, by 

Ue = (^{gevUeiv))'^ 1 + 1 {gevUe{v))' (Jij . (74) 

The group element Uej is fixed, up to a U(l) subgroup, by 

Uefa^n-j = Nljai , (75) 

where for f = {ee') 

JfKf = V4{V) {gev > Ue{v) A Ue'{v))+ * . (76) 

Conversely, every non-degenerate co-tetrad E and spin connection Q give rise to a solution 
via the formulas (71^, (7^, (7^ and f75\ ). 

Proof Let us first consider two consecutive edges e and e' such that / = (ee') . The simplicity 
condition {•kXef{v))-Ue{v) = implies that there exists a 4-vector Nef{v) such that Xef{v) = 
Ue{v)ANf,f{v). Similarly there exists another vector N^'f such that Xe'f{v) = Ue'{v)ANe'f{v). 
The interior closure constraint Xef{v) = Xe'f{v) = Xf{v) requires that Ue'{v) belongs to 
the plane spanned by Ue{v) and Nef{v), so there exist coefficients aef,bef such that 

Ue'{v) = aefNef{v) + hefUe{v) . (77) 
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If Qef = 0, this means that Ue{v) = Ue'{v), since U are normahzed vectors. This is excluded 
by our condition of non-degeneracy, since one would have Xef^iv) AX^'f^^v) = if /2 = (662) 

and = (e'e'a). Denoting Oee' = «7/ therefore has N^j = a^e'Ve' — c^ee'befUe and hence 

X,,,iv) = aee'(t^) (Ueiv) A Ue'iv)') . (78) 

It follows from this expression and the non-degeneracy condition (152|) that the vectors 
Ue{v), e D V, span a 4-dimensional space. As shown in appendix IB 11 the exterior closure 
constraints 

J2^eAv)Xf{v)=0, (79) 

imply the factorization aee' = ^{v)ae{v)ae'{v), where e(u) = ±1 and ae{v) are real numbers, 
independent of the orientation of e, such that 

^aeiv)Ueiv) = and jj = al{v)al{v) sm^Oeeiv) . (80) 

The angle 9ee'{v) is defined by cos^ee'('f) = Ue{v) ■ Ue'{v). These conditions only admit a 
solution if there exists a discrete, geometrical 4-simplex (i.e. a set of edge lengths i{v)) such 
that jf and ^ee' are the areas and dihedral angles in this 4-simplex. In this case, |Q;e(f)| is 
uniquely determined by the spins jf and the unit vectors Ue{v). The ae{v) themselves are 
only fixed up to an overall sign, i.e. if a^, e D f , solves (IHOj) . then e^aeiv), e„ = ±1, is a 
solution as well. 

Given the ordering ei, . . . , 65 of edges at v, we then define 

" \ a (v^ 

ae^Ue^iv),. . . ,aesUerXv)) and Ue{v) = —^==Ue{v) . (81) 

a/I V4I 

These vectors have the property that 

J2 Ueiv) = and Xee'iv) = e{v)V,iv) (Ueiv) A Ueiv)) , (82) 

where e{v) = ±1. Thus, the Ue{v) define tetrad vectors for the 4-simplex dual to v, and we 
can use formula fpB]) to specify corresponding co-tetrad vectors E^iy). 

Next we need to analyze the equations that relate neighboring 4-simplices v and v' , 
connected by the edge e = {yv'): 

Qvv'Ueiv') = Ueiy) , g,,, > Xfiy') = Xfiy) . (83) 

The first condition leads to g^'Ueiy') /\Ue{y')\ = ieUe{v)/\Ue{v)\, where ie = 
sgnae{v) sgnae{v') = ±1. By combining this with the second condition we find that for 
every edge i of the tetrahedron dual to e = {vv') (see appendix IB 2\i 

g,,'EfXv') = e^Eeiv) , (84) 

with the sign e^. = eeSgn{Vi{v)V4{v')) = ±1. We see therefore that the vectors Ei{v) satisfy 
the compatibility condition (iii) in the definition of a co-tetrad, and hence they specify a 
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co-tetrad on the entire simplicial complex. Equation (18^ shows furthermore that g^^i is, up 
to the sign eg, equal to the spin connection 0,^^/ associated with the co-tetrad E: 

Qvv' = ^e^vv' ■ (85) 

In appendix IB 2\ we also derive that the signs e{y) in eq. (!82|) are constant, i.e. e{y) = e(f') 
for neighbouring vertices v and v'. 

The aforementioned ambiguity in the factors ae{v) is transported into the co-tetrad and 
tetrad: for a given solution {jf,nef,Ue,gev,hef), the tetrad and co-tetrad are fixed up to 
a reversal of edges in the geometrical 4-simplices: i.e. up to recplacing {Ee{v),Ue{v)) —>■ 
i-EeXv), -Ue{v)) for all iCv*,eDv. 

If we start, conversely, from a co-tetrad E and its tetrad U, it is clear that equations fl82|) 
and fl85|) define bivectors and connections which solve the equations flSTl) . The associated 
spins jf and group variables Ue and n^f follow directly from the definitions fHTl) and fHSl) . 



A. Determination of h 

So far we have determined {j f,Ue,nef, gev) in terms of a co-tetrad E and signs e and 
ee = e"^""=. In order to complete the characterization of the solution, we also need to 
determine hef. This is done in the following 

Proposition VII. 2 For a non-degenerate co-tetrad E and a choice of global sign e and 
edge signs te, the holonomy of ge = Ce^e around a face f with starting point v has the form 



where the hivector Xf{v) is determined by E as in equation ( [77] ) and Xj = Xf/\Xf\. In 
this equation the bivector is treated as an antisymmetric map acting on and we take an 
exponential of this map. 

In the corresponding solution {jf,Ue,nef,gev,hef) of eqns. l[5l\). the group elements hef 
are uniquely determined, up to gauge transformations, by a choice of angles {9ef, O^f) for 
each wedge: these angles are subject to the conditions 

J20ef = e&f, Y.~^ef = r:Y.ne. (87) 

eC/ eC/ eC/ 

where Qj is the deficit angle of the spin connection. The associated wedge holonomies equal 

Gefiv) = e^^f^f'^-^e^^f^^i^-) (88) 

Proof In order to do the analysis it is convenient to change the frame and base all our 
quantities at the center of the face / (see Fig. H]). That is, we define 

Xefif) = hfe > Xef , Gef{f) = HfeGefkef = hfegevgve'he' f ■ (89) 

It follows from the equations (ISTl) that Xef{f) = Xf{f) is independent of the wedge and 
that Gefif) > Xf{f) = Xf{f). The latter implies that 

Ge/(/) = e'^f^f^^'^e'^f*^f^f^ , (90) 
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Figure 4: The wedge holonomies G(.^f{f) have then' starting and end point at the center of the 
face. 



where Xf{f) = Xf{f)/\Xf{f)\. The angles O^f and O^j have to satisfy constraints, as we 
will show now. 

First we remark that the holonomy around the face / can be written as a product of 
wedge holonomies 

Gfif) = • ■ ■ = ihfe.ge,.) Gfiv) (/l/e. ^7e..) • (91) 

On the right-hand side Gf{v) is the face holonomy based at the vertex v. We have seen 
that the connection of a solution satisfies = egQe, where fl is the spin connection and eg 
is an arbitrary sign. The defining property of the spin connection is that Q^v'Eeiv') = Ei[v) 
for all edges d in the tetrahedron dual to e = {vv'). 

As a result, the holonomy around the face / preserves the co-tetrads associated with the 
triangle dual to /. More precisely, let us suppose that df* = ii + £2 + ^3- Then, the on-shell 
holonomy fulfills 

Gf{v)EeXv)= (^^e}jEeXv), ^ = 1,2,3. (92) 

As we have shown earlier, the bivector Xf{v) is on-shell given by 

Xf{v) = e^{E,,{v)AEe,{v)). (93) 
Hence the condition fl92l) can be equivalently expressed by 

Gf{v)= e'^f^f e"(^^c/ ^<^>^f (^) (94) 

where eg = e™^" and Xf{v) = Xf{v)/\Xf{v)\. The angle 9/ is the deficit angle of the spin 
connection w.r.t. the face /. Combining this result with eq. (!9T|) one obtains that 

eCf eCf ecf 

For a given tetrad U, the associated spin connection Q, and a choice of signs ee, the angles 
9ef and O^f have to meet the constraint fl^ . Once such angles {9ef,9ef) are selected, we 
can solve for hgf recursively. For this, let us set hi = hej and define 

G,(^e./, e,j) = e'^.f^f(-^)e'^J*^i^^^^ . (96) 

The equations 

Ge./(t;0=G.(^e./,4/) (97) 
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can be recursively solved by setting 

hi+l = gei+xViGiQvieihi , hi = hf , (98) 

where hf is an arbitrary initial value. This solution is consistent, since 

hi = hn+l = QeiVnGngvnVn-iGn-l ' ' ' fi'»;2i>i G'l5't,iei ^1 (99) 

deiVnGnidvnVn-lCrn-igVn-lVn) ' ' ' i9v„V„-l ' ' ' 5'l'2 ^1 112 ' ' ' 9vn-lVn)G j ("Wn ) 61 ^1 

= ge^v„G f{Vn)Gj^{Vn)gv„eihi = hi . (100) 

In the third equality, we used that Xef{v') = g^iyXef{v)g~,^. Note that the group element hf 
can be fixed to the identity by a gauge transformation at the face center. This shows that, 
up to gauge, the elements hef are determined by the choice of the angles 9ef, Oej- H 



VIII. SEMICLASSICAL APPROXIMATION OF EFFECTIVE AMPLITUDE 

A. Evaluation of action 

In the previous section, we have seen that solutions of the equations fl48ll51|52l) exist only 
if the set jf is Regge-like and, up to gauge transformation, they are uniquely determined by 
a choice of a discrete metric (coming from a co-tetrad E), of a global sign e, of edge signs 
ee and a choice of U(l) wedge angles {9ef,9ef) subject to flUS]) . 

Proposition VIII. 1 Given a solution characterized by the data {Ue,e,ee,Oef,Of.f), the on- 
shell action is independent of {O^f, O^f) and given by 

Ue^i'^ 7 = 0, 

where Af = (7^ + l^)jf is the area of f in Planck units {SttHG = 1) for 7 > (see eq. 
fTS\) ). Qf is the deficit angle of the spin connection, and Je = (7^^ — 7^) J2fDejf- 

Before giving the proof a few remarks are in order. Firstly, in the EPR model the on-shell 
evaluation is trivial, in agreement with the claim that the EPR model is a quantization of 
the topological sector. Moreover, for general 7, the dependence on the Immirzi parameter 
drops out from the on-shell action. 

Secondly, when evaluating the semiclassical asymptotics of the effective amplitude 
W2(j/)) one has to sum over all classical configurations and hence over eg. This sum gives 
zero unless Je is an even integer. It is interesting to note that when 7^ are both odd integers 
the same condition arises in the spin foam model. 

To see this, note that if 7^ are both odd, the condition that the weight projects down 
to a function of S0(4) (i.e. (7"*" — 7~)j/ G Z) is satisfied without any restriction on jf, 
since (7^ — 7 ) is even. Moreover, the amplitudes in the spin foam model require that the 

invariant SU(2) subspace Inv (^®/DeV^± j is non-trivial. This is the case if and only if J2fjf 

is integer- valued. Therefore, J2fjf integer-valued and Jg is even. 
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Proof As shown in the previous section, the wedge holonomy has the form 
where Xf{v) = Xf{v) /\Xf(v)\. In the SU(2) xSU(2) notation this condition reads 

Recall also that the bivectors Xj^ and X^ are related by 

Xf=j^Xf, \Xf\=jf. (104) 
We insert this into the action, observing that Xj^ /\Xj^\ = 7^/|7^|X^, and obtain 



S = 5^<^2|7+|j/J]ln|tr 



Mi + ^"""^ \g- 

2 I'+M^i 



+ 2|7"|j/$^ln(tr 



/ I ec/ eCf ) 

f e \/De / 



(105) 



(106) 



(107) 



108) 



B. Asymptotic approximation 

In order to arrive at our final result we need to determine the asymptotic approximation 
of the effective amplitude (fl^ for large spins. As shown, this partition function can be 
expressed as an integral 

In = J dx e-^^(") (109) 

over a set of compact variables x. In our case the variables are group elements, so S can be 
taken to be a periodic function. It is customary to restrict the study of this type of integral 
to the case, where S is pure imaginary and use the stationary phase approximation. It is 
less well-know, but nevertheless true, that the stationary phase method is valid when S* is a 
complex function, provided Re(S') > (see [39i] Chapter 7.7). In this reference, it is shown 
that when S is C°° and if |S"p + Re(S') is always strictly positive (with |S"p = d^Sd'^S), 
then the integral is exponentially small. More precisely, if S is C'^''"^ there exists a constant 
C such that 

In < ^—^^:-J—r-^ ■ (HO) 



iV^ mm{\S'\^ + Re{S)) 
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This shows that the integral is exponentially suppressed as long as S" 7^ or Re(5') > 0. 
Therefore, the dominant contribution comes from configurations that are both stationary 



points of the action S, and absolute minima of its real part [39j. One says that is a 
generalized critical point if |S"p(a:c) + Re(S')(xc) = 0. In case there are such points, we have 
the asymptotic approximation 



where Xc are the stationary points of 5*, r is the rank of the Hessian H = didjS{xc), H' 
is its invertible restriction on ker(if)-'- and a is the signature of H'. When the stationary 
points are not isolated, one has an integration over a submanifold of stationary points whose 
dimension equals the dimension of the kernel ker(if). Note that for a generalized critical 
point the action S{xc) is purely imaginary. 

In our case, we have shown that the effective amplitude has no generalized critical points if 
jf is not Regge-like. Then, the previous theorem implies that the amplitude is exponentially 
suppressed. When jf is Regge-like, there is, up to gauge-transformations, a discrete set of 
solutions labelled by {E{jf), eg, e). This result is only valid if one restricts the integration to 
non-degenerate configurations |X A X| > a, with a an arbitrary small positive number. 

When applied to the integral (HM . this gives us that 

W^-^Njf) ~ ^5^exp(iV5l(£;(j,),6,6e)) =^(exp(iV^^) + c.c.), 



(112) 

if the set jf is Regge-like and all Jg even. Otherwise the amplitude is exponentially sup- 
pressed. If there are several tetrad fields E{jf) that correspond to a given set (j/)/ one 
should also sum over them. 

While we have not computed the Hessian, our analysis can give us explicit information 
about its rank r^. In our case, the space of integration is the space of {ue,nef,gev,^ef) 
which is of dimension D = 3E + 2W + 6 x 2E + 6W. Here, E, W, F and V denote the 
number of edges, wedges, faces and vertices of A*. As we have seen, the space of solutions is 
labelled by gauge transformations (Ae,A/,A^) and two U(l) angles {9ef,9ef) subject to one 
constraint per face. Thus, the dimension of the kernel of if is c? = 6E + 6-F + 6V + 2W — 2F. 
We can then compute the rank to be 

rA = D~d = 33E-6V-4:F, (113) 

using the fact that W = AE = lOV. 



C. Degenerate sector 

In order to complete our analysis of the effective amplitude and show its asymptotic 
Regge-like behavior, we have restricted the summation to non-degenerate configurations. 

One could wonder wether the degenerate contributions are dominant or subdominant in 
this semiclassical limit^^. This amounts to asking which sector has the most degenerate 



For instance, in the analysis of the lOj-symbol it was shown that the degenerate configurations were 
non-oscillatory, but dominant 3J, |41 1 
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Hessian, since the amplitude is suppressed by 1/A^2 to the power of the rank of the Hessian. 
Thus, it is the sector with the higher-dimensional space of solutions (higher dimensional 
phase space) that dominates, or in other words the one with higher entropy. 

In order to get an idea of the dimension of the space of solutions in both sectors, let 
us look at the solution of the simplicity and closure constraints at a single vertex. In the 
non-degenerate sector, it is given by 



This describes 10 rotationally invariant degree of freedom, counting 5x4 [/'s subject to 4 
independent constraints minus 6 rotations. These 10 degrees of freedom match the 10 area 
spins. 

On the other hand of the spectrum we can look at the most degenerate contribution, 
where all the wedge products of X's are zero. In this case, the most degenerate solution is 
given by 



Due to the last equation, the Nij are, in effect, 3-dimensional vectors. Now, the counting of 
rotationally invariant degrees of freedom gives 15, 5 more degrees of freedom per vertex than 
in the non degenerate case. Indeed, we have 3 U's plus 3 x 10 A^'s minus 4x3 independent 
constraints minus 6 rotations. 

For each i we can reconstruct a geometrical tetrahedron from Nij,j ^ i, for which N^j 
are the area normal vectors. Hence the degenerate solution determines 5 tetrahedra. These 
5 tetrahedra are "glued together" in the sense that the faces shared by tetrahedra have the 
same area. However, they do not form a 4-simplex. In a 4-simplex the volume of each 
tetrahedron is fixed by the area of the faces, while in the degenerate case the 5 3-volumes 
are independent variables and thus increase the phase space dimension. 

This argument indicates that the phase space dimension of the degenerate configurations 
is higher than the non-degenerate one by at most 5 times the number of vertices. This 
result bears some similarity with the recent canonical analysis of [40], where it was pointed 
out that the phase space dimension associated with spin networks is higher than the corre- 
sponding dimension for discrete geometries. Our reasoning suggests that this extra phase 
space corresponds to 4d degenerate solutions. 

This analysis is suggestive, but not complete, since one would need to analyze the gluing 
equations and the other degenerate sectors. However, it leads one to suspect that the de- 
generate configuration dominate the effective amplitude in the semiclassical limit if they are 
included. In this case, the non-degeneracy requirement would be necessary. One challenge 
is to be able to formulate this requirement at the level of the spin foam model and not only 
in the path integral representation. Another possibility is that the degenerate contributions 
are suppressed when we couple the effective amplitude to a semiclassical boundary state. 
This is a scenario that has been realized in the case of the lOj-symbol jl]. 




(114) 





Ni.-U = 0. 



(115) 



IX. SUMMARY AND DISCUSSION 



In this work, we have studied the semiclassical properties of the Riemannian spin foam 
models FK7. We have shown that, in the semi-classical limit, where all the bulk spins are 
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rescaled, the amplitude converges rapidly towards the exponential of i times the Regge ac- 
tion, provided the face's spins can be understood as coming from a discrete geometry. When 
the spins do not arise from a discrete geometry, the spin foam amplitude is exponentially 
suppressed. 

There are several remarks to be made about this result: First, it is shown for an arbitrary 
triangulation and not only for the amplitude associated with a 4-simplex. This should be 
contrasted with what was achieved in the context of the Barret-Crane model, where only 
one or two 4-simplices were considered |8||. An extension of these results to more 4-simplices 
seemed increasingly complicated (see [4j] for a very recent discussion of this in the context of 
Regge calculus). The second fact to be noticed is that the Immirzi dependence drops out in 
the semiclassical limit. This should indeed be the case, since nothing depends on the Immirzi 
parameter at the classical level (except when it is zero). Nevertheless it was not obvious 
from the original definition of the amplitude that this would happen. Also, the results shown 
here depend heavily on the details of the implementation of the simplicity constraints: they 
rely on the specific choice of the measure Z)J^ (see eq. ([3])). For instance, we cannot extend 
our results to the ELPR7 model for 7 > 1 which includes the Barret-Crane model for 
7 = 00. A fourth point concerns the fact that in spin foam models areas are the natural 
variables, whereas one needs access to edge lengths in order to have a discrete geometry. To 
formulate constraints on areas, so that they correspond to discrete geometries, has been so 
far one of the conundrums faced in the LQG/spin foam ap pro ach. Several studies have been 



launched in order to tackle this problem (see for instance [331, |3J, l35|) , but the results show 
that performing this explicitly is an incredibly difficult algebraic task. What we find quite 
remarkable is that it is not necessary to answer this question analytically to get the proper 
semiclassical limit of a spin foam model. The spin foam model "knows" which set of areas 
does or does not arise from a 4d geometry and it naturally suppresses the non-geometric 
phase in the semiclassical limit. 

These results provide considerable evidence in favor of the proposed spin foam amplitude 
as a valid amplitude for quantum gravity, in the sense that it reproduces expected semiclas- 
sical behavior. There is, however, more work to be done to fully confirm this picture. 

First of all, in order to obtain this result we have to restrict the summation to non- 
degenerate configurations. We know how to implement this restriction in the path integral 
formulation, but not in terms of the spin foam model. As we have argued, this restriction 
may be important in order to get the correct semi-classical limit, but a deeper analysis is 
clearly required to establish this firmly. 

More crucially, we have shown the semiclassical property of the bulk amplitude, where the 
bulk spins are fixed and uniformly rescaled to large values. That is, we have demonstrated 
the proper semiclassicality for certain histories that one should sum over in computing 
amplitudes. What we are ultimately interested in is the semiclassical property of the sum 
over amplitudes. Given a boundary spin network, we would like to sum over all spins in the 
interior compatible with the boundary spin network and show that the resulting amplitude 
gives an object that can be interpreted as the exponential of the Hamilton-Jacobi functional 
of a gravity action. Our result is a necessary condition for this to happen, but we have not 
shown that this is sufficient. 

What would be required is that for given semiclassical boundary states peaked on large 
spins, the corresponding amphtude is peaked around large bulk spins as well; and that the 
semiclassical amplitude reduces effectively to a summation over discrete geometries with the 
Regge action. In a sense, one needs that the large spin limit and the integration over the 
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spins commute with each other. Whether this happens or not is not obvious: one might be 
worried, for instance, that the summation over spins is much less restricted than a summation 
over discrete geometries and that this will lead to stronger equations of motions. It might 
be, on the other hand, that the exponential suppression of non Regge-like configurations 
is strong enough to effectively reduce the summation to a sum over geometries. This is an 
important question that deserves to be studied further. 

An obvious open problem is whether our results can be extended to the Lorentzian case. 
We expect that this is possible, however, it has not been shown yet wether the present 
Lorentzian models admit a nice action representation, which is needed for our analysis. 

Moreover, our work does not address the question of the continuum limit of spin foam 
models. We have considered the semiclassical limit of discrete configurations on a fixed 
triangulation. One might want to take a continuum limit, where the number of bound- 
ary vertices of the spin network grows. It is not clear if such a limit commutes with the 
semiclassical limit taken here. 

Despite all these open questions, we feel that the semiclassicality shown here opens the 
way towards new, exciting developments in the spin foam approach to quantum gravity. 
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Appendix A: RELATION BETWEEN CO TETRAD AND TETRAD IN A 

4 SIMPLEX 

Based on the duality (!63|) between discrete tetrad and co-tetrad, we can prove a number 
of identities that are analogous to equations for the co-tetrad and tetrad in the continuum. 
Consider a vertex v in A* and label the vertices p G v* (and corresponding dual edges e D v) 
by lowercase letters i,j,k... = 1, . . . , 5. The tetrad and co-tetrad vectors Ue{v) and Eee'{v) 
are written as f/* and Eij. We denote M'^-indices by capital letters J, J, K etc. 

Proposition A.l The equation 

U}E^,' = 51-61 (Al) 

determines a bijection between non-degenerate vectors Eij G M.^, i, j = 1, . . . ,5, satisfying 
E,j+Eji = 0, Eij + Ejk + Eki = Vz,j,fc = 1,...,5, (A2) 

and non-degenerate vectors W E , i = 1, . . . ,5, for which Yl^=i f^* = 0. 
The map from Eij to W is given by 

= ^ E ^ (^nk A E,,, A E,3,) , (A3) 
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where k is any vertex different from i. W is independent of this choice thanks to the identity 
M^) . V4 denotes the oriented volume of the 4^parallelotope spanned by the co-tetrad vectors, 

V^ = det{E2i,...,E,,) , (A4) 

and we have set 

[^{E, A • • ■ A i?„)]/,../,_„ = ei,...jX'-" ---Et. (A5) 
The norm of W is proportional to the volume V3 of the tetrahedron orthogonal to W : 

\U'\ = ^ (A6) 

The inverse of V4 equals the determinant of the U 's: 

-^ = det([/2i,...,f/5i) (A7) 

The inverse map from U to E is specified by 

More generally, the relation between U and E is given by 

(A9) 

r/ie special cases n = 1 and n = 3 return equations ^A3\} and ^A8\} respectively. For n = 2 
one obtains 

V,U\jWj] = 5^e^*^"^'^^e^V7vi?,n/Xfc'^, k ^ t,j . (AlO) 



Proof For the first part of the proof, let us assume the vectors Ei^ with property (lA2p 
are given and that the t/j's satisfy relation (lAip . The identity (1A3P is proven, like in the 
continuum, by contracting the left- and right-hand side with Ej^. That the t/j's close follows 
directly from (lA3p . Formula (1A6P can be derived by using (1A3P and the relation between 
volume and Gram's determinant. Identity flA7p follows from flAip and the multiplication 
rule for determinants. By contraction and use of flA7p . we also verify eq. flASp . 

To demonstrate that the right-hand side is independent of k i, it helps to regard the 
Eik as edge vectors of a 4-simplex in W^. We can think of this 4-simplex as the image of the 
4-simplex a C A under an affine transformation. Let Pi, . . . , P5 G denote the images of 
the vertices pi, . . . ,^5 C cr. Then, the edge vectors are equal to 

E,,fc = P,-Pfc. (All) 

Without loss of generality, we can assume that 

5 



1=1 
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Using this, we deduce that 

U' = ^^ Yl e'''-'^''' HPn ^ Pn ^ Pn) . (A13) 

making the independence oi k ^ i in (lASp manifest. 

Conversely, suppose we have vectors Ui that close and that the -Ejfc's fulfill relation (lAip . 
We then define vectors 

^. = ^^4 ^kJ^,^.^^^{U'' ^U'' AU'') (Al4) 

and verify that 

E,k = Pi-Pk- (A15) 

Hence the vectors Eik close. 

Relation (lA9p is demonstrated by contracting with n E's. M 

Appendix B: RECONSTRUCTION OF 4 GEOMETRY 

In this appendix, we complete the proof of proposition fIVII.ll) . In the first part, we will 
derive that the bivectors Xf{v) arise from a geometric 4-simplex. A key step for this is that 
the factors Oee' in eq. ( ITHl) factorize. In the second part, we derive relations among tetrad 
vectors between neighbouring vertices, showing that the tetrad and co-tetrad vectors define 
a consistent discrete geometry on the simplicial complex. We will prove, in particular, that 
the sign factors e{v) in eq. fl5^ are the same for every vertex. 

1. Reconstruction of 4 simplex 

Consider a vertex t> C A* and the edges ei, . . . , 65 D f . To simplify formulas, we use the 
abbreviations Xij = X^^^j , Ui = U^^ and Eij = E^^e^ ■ 

Proposition B.l Let Xij = —Xji, i,j = 1, ... ,5, be non-degenerate bivectors (i.e. \Xij A 
Xki\ > 0) which satisfy the simplicity and closure constraint 

x^/mj = o, (Bi) 

^X,,=0. (B2) 

Then, there are, modulo translations, precisely two 4^simplices whose area bivectors equal 
•kXij and they are related by a reversal of edge vectors. That is, there are exactly two sets of 
vectors Eij G M^, i, j = 1, . . . ,5, obeying the closure condition M^j) . such that 

■^"■i ^ ^ 2 {P'nik A Enk) , k ^ i,j . (B3) 

m,n 

The sign factor e is either 1 or —1 W i,j = 1, . . .5. The two sets {Eij} are related by the 
30(4.) transformation Eij —Eij- 
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Proof The simplicity constraints (IBip imply that 

X,j = a,jUi A Uj , (B4) 

where Gl i^j IS 3j symmetric matrix of normalization factors and the wedge product stands for 
the bivector 

{U, A = Ui^'U/^ = U/U/ - il/U,' . (B5) 

The closure constraint states that 

^aijlJi^Uj = U,^^aifij = Q Vi = l,...,5. (B6) 

Consequently, 

5 

i=i 

for suitable diagonal elements an. 

Next we eliminate one of the five Uj in the last equation, say, JJm- For arbitrary k, /, 

{akmOtij - aimOtki) Uj = E {otkmaij - aimOikj) Uj = 0. (B8) 

Since the bivectors are non-degenerate, four of the five normal vectors f/j must be linearly 
independent. Therefore, 

akmOij = akjOim ■ (B9) 

In particular, for I = j, 

akmOijj = akjOjm . (BIO) 
By non-degeneracy, all a^j are non-zero, so 



ttfcm = = — . (Bll) 



Let us pick one j = jo and define 



a = -^M= . (B12) 



a 
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Then, 

= sgn{ajgjg) aiaj (B13) 

and the bivectors have the form 



Xij = I [oiiUi^ A [oLjUj^ , (B14) 
where e = sgn(Q;j(,j(,) is a sign independent of i and j. From eq. (]B7p we also know that 



3 

By taking the square of eq. (lB14p . we get 



3ij — Q^i ctj sin^ , cos 9ij = U ■ Uj , (B16) 
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Figure 5: Choice of labelling at neighbouring vertices v and v' . 



which fixes the modulus of given jij and f/j. Eq. ( IBISI) implies furthermore that the signs 
sgnoj are fixed up to an overall sign change ai —Qi, i = 1, . . . , 5. 

At this point, we can reconstruct the tetrad and co-tetrad vectors. Define 



OLiUi 



with V4 = det ( a2t^2, • • • , ot^U^ ) • 



Then, we obtain that 



and 



- = det(?72,...,t/5) 



(B17) 



(B18) 



Xij = e II/4I f/i A Uj = tViU^A Uj , (B19) 
where e = esgn(V4). By proposition I VI .41 and I A . 11 the f/j's define corresponding dual vectors 

^ ^ 9 ^kijmn * (-E^mfc A Enk) , k i,j . (B20) 



Eij such that 



2. Reconstruction of co tetrad and tetrad 

Next we deal with the equations fl83p that relate variables from neighbouring 4-simplices. 
We consider an edge e = {vv') and employ the following shorthand notation: 

Uo = U,{v), U'^ = g,,,U,{v'), (B21) 

U, = U,^{v), Ui = g,,,U,'^{v') , (B22) 

E,, = (v) , E',j = w^e.e, {v') . (B23) 

The labels i are chosen such that (cjee^) corresponds to one of the four faces adjacent to 
e (see fig. [5]). One can check that this ordering is compatible with our requirement that 
orientations of neighbouring 4-simplices are consistent. 

As seen in section IVIH the exterior closure constraints leads to 



Y^U, = -Uo and ^f/; = -?7^ 

i i 



(B24) 
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Moreover, due to the eqns. (183|) . the U and U' are related as follows: 

Xoi = eViUo A a,) = eV'(t/^ A U^) , (B25) 



K\ Vol' 



where e, e', e = ±1 and 

1/V = det (Ui, U2, Us, Ui) , l/V = det iU[, U^, U'^, U'^) . (B26) 

Proposition B.2 The conditions \B24\ \B25\) imply that 

e = e', i=sgn{VV')a, VUq = aVU^^ and U- = aU, + a,Uo , (B27) 

where a is an arbitrary sign factor and ai are coefficients such that Yli'^i = « (l — y?)- 
Moreover, for the co-tetrad vectors Eij and E[- one has the identity 

Elj = aEi, . (B28) 

Proof The equations (1B25I1 tell us that [/q is proportional to Uq and that f// is a linear 
combination of Ui and Uq. More precisely, 

t/; = iee'l^^U^ + aiUo , (B29) 

Using the relation (lB29p . we obtain 

1/V" = det ([/;, (7^, (7^, U',) = det ([/; (7(, C/^, C/^) (B30) 

Thus, e = e' and |f/o|^ = il^ol^'- By defining the sign factor 
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v 



a ^ ~ei-^ , (B32) 



we arrive at eq. ( 1B27[) . 

By using eq. flASp of prop. flA.ip . we can now compute explicitly the relation between 



co-tetrad vectors for edges that are shared by the 4-simplices dual to v and v': 
1 
3! 



^jk = Iv'e.k'''''' * ^ K ^ K) = c^'Ve.k'''''' * (f/n A U,, A U,,) = aE,k ■ (B33) 



The relation E^j = aEij shows that the Eij satisfy the metricity condition (iii) in the 
definition of a co-tetrad. Therefore, the co-tetrad and tetrad vectors determine a consistent 
4-geometry on the simplicial complex. 
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